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Abstract 

We apply the topology conserving gauge action proposed by Liischer to the four-dimensional 
lattice QCD simulation in the quenched approximation. With this gauge action the topological 
charge is stabilized along the hybrid Monte Carlo updates compared to the standard Wilson gauge 
action. The quark potential and renormalized coupling constant are in good agreement with 
the results obtained with the Wilson gauge action. We also investigate the low-lying eigenvalue 
distribution of the hermitian Wilson-Dirac operator, which is relevant for the construction of the 
overlap-Dirac operator. 

PACS numbers: 
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I. INTRODUCTION 



Chiral symmetry and topology are tightly related with each other in the gauge field 
theory through the quantum correction. Namely, the axial anomaly appears at the one-loop 
level, and its integral over space-time leads to the topological charge of the background 
gauge field. In principle, one should be able to analyze the implication of this relation for 
physical observables, such as the neutron electric dipole moment, using the lattice gauge 
theory, which provides a rigorous formulation of the non-Abelian gauge theories even in the 
non-perturbative regime. Such study is very difficult with the Wilson-type Dirac operator, 
since the chiral symmetry is explicitly violated on the lattice. 

The overlap-Dirac operator |l[ 

D = -[l + -f5Sgn{aHw)] , a = — (1) 
a I + s 

realizes the exact chi.a, =y..et^ at fl,.te lattice spacing . B satisfying the Ginspa.g- 
Wilson relation |^ 

-f,D + D-f, = aD-f,D. (2) 

It is constructed from the Wilson-Dirac operator aDw with the Wilson parameter r = 1; 
the hermitian Wilson-Dirac operator aHw = 75(a-Dvi/ — 1 — s) enters as an argument of the 
sign function sgn(x). The parameter s in is a fixed number in the region \s\ < 1. 

Since the definition contains a non-smooth function, the locality of the Dirac operator 
could be lost when there are near-zero eigenvalues of |aifvK|. This is consistent with the index 
theorem, because the index of the Dirac operator, which may be considered as a definition 
of the topological charge, is a non-smooth function of the background gauge field. When 
the topological charge changes the value, the Dirac operator must become ill-defined, and 
this is exactly the point where aHw has a zero eigenvalue. 

The locality of the overlap-Dirac operator is guaranteed for the gauge fields on which 
;he minimum eigenvalue of |aifvy| is bounded from below by a positive (non-zero) constant 
^. This condition is proved to be satisfied if the gauge field configuration is smooth and 
each plaquette is close enough to one; 



||1 - P^^(x)|| < e for all x, (/i, z/). (3) 
Here, P^u{x) is the plaquette variable at x on the fi-u plane, and || ■ ■ ■ || denotes the norm 
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of the operator. In the four-dimensional case, the parameter e ~ 1/20.49 is a sufficient (but 
not a necessary) condition Qj. This is called the "admissibility" bound. 

One can construct a gauge action, which generates gau^e configurations respecting the 
condition (jSj). For instance, Liischer proposed the action 

( l-ReTrP^.(x)/3 

V p p ^ ^/Q^/ ' ^^^^ 1 ~ ReTrP^^ X /3 < e, 

= <^ p^-i^- ReTrP^^(a;)/3)/e ^ (4) 

[ 00 otherwise 

which has the same continuum limit as the standard Wilson gauge action does. In fact, 
the limit e = 00 corresponds to the standard Wilson gauge action. Unfortunately, the 
bound e ~ 1/20.49 is too tight to produce gauge field ensembles corresponding to the lattice 
spacing around 0.1 fm; for practical purposes, one must choose much larger values of e. 
An interesting question is, then, whether the action can keep the good properties for e 
significantly larger than 1/20.49. To be explicit, one expects that (i) the topology change 
during the molecular-dynamics-type simulation is suppressed, and (ii) the appearance of the 
near-zero eigenvalue of |a-f^iy| is suppressed, compared to the standard Wilson gauge action. 
The point (i) is important in order to efficiently generate gauge configurations with large 
topological charge, which is necessary for the study of the e-regime or the ^-vacuum. With 
the point (ii), the locality of the overlap-Dirac operator is improved, and the numerical cost 
to apply the overlap-Dirac operator is reduced. In the numerical application to the massive 
Schwinger model with e = 1, the stability of the topological charge and the improvement 



of the chiral symmetry with the domain-wall fermion were observed 



Also, in the 



bur- dimensional quenched QCD, good stability of the topological charge has been reported 

For gauge actions to be useful in practical simulations, good scaling property toward the 
continuum limit is required. The action differs from the standard Wilson gauge action 
only at O(a^), and we expect that it approaches to the continuum limit as quickly as the 
standard action does. The scaling would be better for the improved gauge actions, such 
as the Liischer- Weisz Iwasaki or DBW2 [igI] gauge actions, but an advantage of 

is that it contains a parameter which directly controls the admissibility bound and thus 
the appearance of the low-lying modes of aHy/- For the improved actions including the 
rectangle loop, on the other hand, the low-lying modes are suppressed for large values of the 
rectangle coupling {e.g. with the DBW2 action) at the price of loosing the good scaling for 



short distance quantities [17 1. 

The goal of this paper is to give a systematic quenched QCD study of the topology 
conserving gauge action. We find that the topology change is indeed suppressed when the 
parameter e is of order one. We also find that the scaling violation in the static quark 
potential remain reasonably small and the tadpole improved perturbation theory for the 
renormalized gauge coupling shows a good convergence in the parameter range of our study. 
Therefore, the topology conserving gauge action has a desired properties for a practical 
application. Here we would like to comment on the possible application for the future work 
with the dynamical overlap fermion. In the standard method, with the Wilson plaquette 
action, one projects out the smallest eigenmodes of at every molecular dynamics step in 
the simulation trajectory and judge if the topology change occurs or not. When the topology 
change occurs, one recalculates the link updation with much higher accuracy on the topology 
crossing point to choose either entering the new sector (refraction) or going back to the 
previous sector (reflection) Q]. If one omits this step it would make the acceptance very 
low due to the non-smoothness of the determinant as a functional of the gauge configuration. 
On the other hand, if one uses the topology conserving action, crossing Hw = can be 
strongly suppressed, so that one can avoid the CPU time consuming reflection/refraction 
method. In this sense, the use of the topology conserving gauge action can also be useful in 
full QCD simulations with the overlap fermions. Combination with the stout link version of 
the overlap fermion [3] would be interesting as well. 

This paper is organized as follows. After describing the simulation methods in Section ITU 
we show the fundamental scaling studies in Section lTTH that is the determination of the lattice 
spacing and a scaling test with the static quark potential. Renormalization of the coupling 
constant with the action can be estimated using perturbation theory as described in 
Section IIVI Section |3 is the main part of this paper; we report how much the topology 
change may occur with different choices of parameters. In Section IVII the locality and 
the numerical costs of the overlap fermion with gauge fields satisfying the bound are 
discussed. Conclusion and outlook are given in Section I VI II 
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II. LATTICE SIMULATIONS 



Although several types of the gauge action that generate the "admissible" gauge fields 



satisfying the bound are proposed [10|, we take the simplest choice (jH). We study 



2l|. 



three values of 1/e: 1, 2/3, and 0. Note that 1/e = corresponds to the conventional Wilson 
gauge action. The value 1/e = 2/3 is the boundary, below this value the gauge links can 
take any value in the gauge group SU{3) and the positivity is guaranteed j20|. 

The link variables are generated with the standard hybrid Monte Carlo algorithm 
We take the molecular dynamics step size At in the range 0.01-0.02 and the number of 
steps in an unit trajectory Nmds = 20-40. During the molecular dynamics steps we monitor 
that the condition 1 — ReTrP^i,(a;)/3 < e is always satisfied with our choice of the step size. 
We discarded at least 2000 trajectories for thermalization before measuring observables. 

In order to measure the topological charge, we develop a new type of cooling method. It 
consists of the hybrid Monte Carlo simulation with an exponentially increasing (3 value /^cooi 
and decreasing step size Atcooi as a function of trajectory rit, i.e. 
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/3cool = Anit X (1.5)' 

Ar,ooi = Ari,it x (1.5)-"*/^, (5) 

with a fixed 1/ecooi- Note that \/ Aooi AtcooI is fixed so that the evolution at each step is kept 
small. This method allows us to "cool" the configuration smoothly, keeping the admissibility 
bound © with 1/e = 1/ecooi- For the parameters, we take (Anit, Arinit, l/ccooi) = (2.0, 0.01, 
1) for the configurations generated with 1/e = 1, and (3.5, 0.01, 2/3) for the configurations 
with 1/e = 2/3 or 1/e = 0. Even for the gauge configuration generated with the standard 
gauge action (1/e = 0), the condition 1/ecooi = 2/3 can be used because it allows all values 
of SU(3). After 50-200 steps, the link variables are cooled down close to a classical solution 



in each topological sector. In fact, the geometrical definition of the topological charge 
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Qgeo ^i^Y. e^'^'^ReTr (P^.(x)P,.(x)) (6) 

X 

of these "cooled" configurations gives numbers close to an integer times a universal factor 
Z~}^. Namely, Q = ZgeoQgeo is close to an integer. We determine Zgeo through would-be 
Q = I gauge configurations, as Z~J-^ = 0.923(4). As Figure [T] shows, the topological charge Q 
is consistent with the index of the overlap-Dirac operator with s = 0.6, which is calculated 



L=16, beta=1.42, 1/e=1 



L=16, beta=6.13, 1/e=0 



j 


1 1 1 1 1 1 

the overlap index — \ — 
geometorical Q ^ 












2110 2120 2130 2140 2150 2160 2170 2180 2190 2200 
conf 



8000 8100 8200 8300 8400 8500 8600 8700 8800 8900 9000 
conf 



FIG. 1: Comparison of the topological charge Q, calculated from the geometrical definition Qgeo 
(jnj as Q = ZgeoQgeo, with the index of the overlap-Dirac operator. Qgeo is obtained by the cooling 
method (stars), and the index (pluses) is calculated for the overlap-Dirac operator with s = 0.6. 
The agreement is better for 1/e = 1 (left) than for 1/e = (right). 

as described in Section IVII The consistency is better for 1/e = 1 than for the standard 
Wilson gauge action 1/e = 0. 

To generate topologically non-trivial gauge configurations, we start the hybrid Monte 
Carlo simulation with the initial condition 

/ ^2TTiX4Q/L'^ \ ( 

U2{X) = 



\ 

(l 



\ 



-2iTix2/L'^ 



J 



\ 



g27rja;3(52,2,L-l/-f' 



^-2TTixiQ5^^^L_i/L 



\ 



(7) 



which is a discretized version of the classical solution on a four-dimensional torus 2^. It 
can be used for any integer value of Q. We confirmed that the topological charge assigned 
in this way agrees with the index of the overlap operator with s = 0.6. 

We summarize the simulation parameters and the plaquette expectation values (for the 
run with the initial configuration with Q = 0) in Table IH The length of unit trajectory 
is 0.2-0.4, and the step size is chosen such that the acceptance rate becomes larger than 
~ 70%. 
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Lattice size 1/e (3 At Nmds acceptance plaquette 



1 


1.0 


0.01 


40 




1.2 


0.01 


40 




1.3 


0.01 


40 


2/3 


2.25 


0.01 


40 




2.4 


0.01 


40 




2.55 


0.01 


40 





5.8 


0.02 


20 




5.9 


0.02 


20 




6.0 


0.02 


20 


1 


1.3 


0.01 


20 




1.42 


0.01 


20 


2/3 


2.55 


0.01 


20 




2.7 


0.01 


20 





6.0 


0.01 


20 




6.13 


0.01 


40 


1 


1.3 


0.01 


20 




1.42 


0.01 


20 


2/3 


2.55 


0.01 


20 




2.7 


0.01 


20 





6.0 


0.015 


20 




6.13 


0.01 


20 



89% 0.539127(9 

90% 0.566429(6; 

90% 0.578405(6; 

93% 0.55102(1) 

93% 0.56861(1) 

93% 0.58435(1) 

69% 0.56763(5) 

69% 0.58190(3) 

68% 0.59364(2) 

82% 0.57840(1) 

82% 0.59167(1) 

88% 0.58428(2) 

87% 0.59862(1) 

89% 0.59382(5) 

88% 0.60711(4) 

72% 0.57847(9) 

74% 0.59165(1) 

82% 0.58438(2) 

82% 0.59865(1) 

53% 0.59382(4) 

83% 0.60716(3) 



TABLE I: Simulation parameters and the plaquette expectation values (for the run with the initial 
configuration with Q = 0). 

III. STATIC QUARK POTENTIAL 

In this section we describe the measurement of the static quark potential to determine 
the lattice spacing for each parameter choice. We then compare the scaling violation and 
the rotational symmetry violation with the case of the standard Wilson gauge action. In the 
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following, we assume that the topology of the gauge field does not affect the Wilson loops, 
and choose the run with Q = initial configuration for the measurement. 

We measure the Wilson loops W{f,t) using the smearing technique according to 
where the spatial separation f/a is taken to be an integer multiples of elementary vectors 
V = (1,0,0), (1,1,0), (2,1,0), (1,1,1), (2,1,1), (2,2,1). With the assumption that the 
Wilson loop is an exponential function for large temporal side t/a, (W^(r, t)) = exp{—V{r)t), 
we extract the static quark potential aV{r). The measurements are done every 20 trajectories 
and the errors are estimated by the jackknife method. 

As a reference scale, we measure the Sommer scales ro and Vc defined as ro-F(ro) = 

1.65 and rlF{r^ = 0.65, respectively. Here, the force F{r) on the lattice is given by a 
derivative in the direction of u/a = (1, 0, 0); 

\u/a\ 

for f/\f\ = (1, 0, 0). r/ is introduced to cancel the discretization error in the short distances, 
using the one-gluon exchange potential on the lattice 

1 aG{f) — aG{f — u) 

4:7!-{rj/ay \u/a\ ' 

r 'Pk n-.,<!os(oAj/a) 
'"^"•'-l(2^ 4E,t.Wfe/2) - 

In Table HTl we list the values of the Sommer scales ro/a, rc/a as well as their ratio Tc/tq. 
The numerical results for aV{f) and r']F{ri) for the case that f/a is an integer multiples of 
u/a are given in Tables UTTl and HVl The values of rj/a are also listed. 

The scaling can be tested for the ratio Vc/tq. Figure |21 presents the dependence of this 
ratio for different values of 1/e. Our results for 1/e = 2/3 and 1 are in perfect agreement 
with the pre vious high statistics study for the standard Wilson gauge action by Necco and 
Sommer [2^. Moreover, we do not find any statistically significant scaling violation except 
for the coarsest lattice points around (a/ro)^ ~ 0.1. 

Figure 01 shows a comparison of the potential itself in a dimensionless combination, i. e. 
^o^(^) = f^oi^i'f^) ~ y{f^c)) versus |r|/ro. For V{rc) we interpolate the data in the direction 
f/r= (1, 0, 0). The data at (3 =1.3, 1/e = 1 are plotted together with the curve representing 
the continuum limit obtained in [23|. The agreement is satisfactory (less than two sigma) 
for long distances r/ro > 0.5. 



Lattice size 


1/e B 


statistics 


Tn/ci 
' u/ ^ 


TrJa Tn/vn 

' C/ ^ ' C/ ' \J 




1 1.0 


3800 


3 257('30) 


1 7081f50) 5244('52) 




1.2 


3800 


4 555(73) 


2 319fl0) 5091(81") 




1.3 


3800 


5.140(50) 


2.710(14) 0.5272(53) 




2/3 2.25 


3800 


3.498(24) 


1.8304(60) 0.5233(41) 




2.4 


3800 


4.386(53) 


2.254(10) 0.5141(61) 




2.55 


3800 


5.433(72) 


2.809(18) 0.5170(67) 


16^ 


1 1.3 


2300 


5.240(96) 


2.686(13) 0.5126(98) 




1.42 


2247 


6 240(RQ) 


3 970l'26'l 5941 ('83'! 




2/3 2.55 


1950 


5.290(69) 


2.738(15) 0.5174(72) 




2.7 


2150 


6.559(76) 


3.382(22) 0.5156(65) 


continuum limit [28] 






0.5133(24) 



TABLE II: Sommer scales rg/a, Tc/a and their ratio. 

For short distances, on the other hand, we can see deviations of order 10%, as shown 
in Figure m where a ratio {y{r) — V'cont(|^))/Vcont(|^) is plotted. Kontd'l) represents the 
curve in the continuum hmit drawn in Figure El The points corresponding to the separation 
r/a = (1,0,0) and (2,0,0) deviates significantly from zero in the upward direction, while 
the points (1, 1, 0) and (1, 1, 1) are lower than zero. This implies the rotational symmetry 
violation. Figure El (left panel) shows the size of the rotational symmetry violation at the 
point (1,0,0) as a function of the lattice spacing. We find that the size of the violation is 
quite similar for different values of 1/e including the standard Wilson gauge action. It does 
not show a tendency that the rotational symmetry violation goes to zero in the continuum 
limit, but this makes sense because the relevant scale of the observable is also diverging as 
1/a. After correcting the tree level violation by introducing dj as l/(47r(i/) = G{d), which 
is an analogue of r/ in but is defined for the potential, we obtain the plot on the right 
panel of Figure El It is indeed improved. The remaining correction is of order ^^(l/a), 
which vanishes as ~ l/ln(l/a) near the continuum limit. 

These observations are consistent with the fact that the topology conserving gauge action 
has the same O(a^) scaling violation as the standard Wilson gauge action. The difference 
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1/e = 1 

(3 r/a rj/a 


12^ 

aV{r) rjFin) 


16^ 

aV{r) rjF{ri) 


1.0 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 


0.50459(20) 

0.77828(61) 0.5056(10) 
0.9629(15) 0.9520(69) 
1.1176(27) 1.691(26) 
1.2623(45) 2.751(80) 
1.4052(77) 4.33(22) 




1.2 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 


0.44877(16) 

0.65982(39) 0.38993(65) 
0.78291(80) 0.6346(34) 
0.8775(13) 1.034(10) 
0.9588(29) 1.545(45) 
1.0322(47) 2.23(12) 




1.3 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 


0.42730(10) 

0.61711(34) 0.35252(99) 
0.72140(69) 0.53909(48) 
0.7977(12) 0.848(14) 
0.8608(21) 1.240(36) 
0.9230(25) 1.887(85) 


0.42709(20) 

0.61710(66) 0.35099(68) 
0.72130(92) 0.5490(29) 
0.7961(15) 0.8325(81) 
0.8583(23) 1.180(32) 
0.9150(27) 1.809(79) 
0.9636(51) 1.93(24) 
1.0215(51) 3.09(37) 


1.42 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 




0.40443(15) 

0.57416(43) 0.31444(58) 
0.66091(75) 0.4567(22) 
0.7200(12) 0.6583(61) 
0.7691(17) 0.940(14) 
0.8076(24) 1.189(48) 
0.8457(30) 1.675(64) 
0.8832(37) 1.91(14) 



TABLE III: Potential and force values for the case that r/a is an integer multiples of the unit 
vector u/a = (1,0,0). Results for 1/e = 1. 
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1/e = 2/3 
f3 r/a rj/a 


12^ 

aV{r) rjF{ri) 


16^ 

aV{f) rjF{rj) 


2.25 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 


0.48470(15) 

0.74012(57) 0.47189(97) 
0.9077(13) 0.8640(56) 
1.0463(22) 1.515(21) 
1.1701(38) 2.353(64) 
1.2901(58) 3.64(15) 




2.4 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 


0.44908(12) 

0.66434(41) 0.39770(70) 
0.79152(84) 0.6557(37) 
0.8889(15) 1.065(12) 
0.9749(23) 1.635(32) 
1.0541(30) 2.401(74) 




2.55 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 


0.42013(11) 

0.60682(36) 0.34493(58) 
0.70826(72) 0.5230(28) 
0.7806(13) 0.7913(90) 
0.8430(18) 1.187(18) 
0.8986(23) 1.686(37) 


0.42042(16) 

0.60786(51) 0.34590(72) 
0.71227(95) 0.5337(32) 
0.7878(16) 0.8211(93) 
0.8538(22) 1.210(21) 
0.9157(29) 1.765(47) 
0.9710(43) 2.229(84) 
1.0266(52) 2.94(15) 


2.7 1.0 

2.0 1.36 
3.0 2.28 
4.0 3.31 
5.0 4.36 
6.0 5.39 
7.0 6.41 
8.0 7.43 




0.39590(15) 

0.56100(44) 0.30650(53) 
0.64733(62) 0.4456(22) 
0.70527(90) 0.6329(56) 
0.7528(14) 0.907(14) 
0.7937(19) 1.309(28) 
0.8321(24) 1.531(44) 
0.8703(29) 2.035(80) 



TABLE IV: Same as Table |lTll but for 1/e = 2/3. 
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FIG. 2: A ratio of the Sommer scales Vc/vq. Squares and triangles are data for the topology 
conserving gauge action with 1/e = 2/3 and 1, respectively. Open circles represent the standard 



Wilson gauge action (from 



281 ]) and the filled circle is their continuum limit. 



starts at 0{a^)^ which is not visible at the level of precision in our numerical study. 

Finally, we confirm our assumption that the topology does not affect the quark potential 
by measuring rg for two initial value of Q (0 and —3). Measurements are done on a 16^ 
lattice aX (3 = 1.42, 1/e = 1, for which the probability of the topology change is extremely 
small as discussed in the next section. Our results are ro/a = 6.24(9) for the (5 = initial 
condition and 6.11(13) for Q = —3. 
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beta=1.3, 1/e=1 



o 

> 

I 

> 

o 




unitvec=(1 ,0,0) 
unitvec=(1 ,1 ,0) 
unitvec=(2,1,0) 
unitvec=(1 ,1 ,1) 
unitvec=(2,1 ,1) 
unitvec=(2,2,1) 
continuum limit (Necco,Sommer 01) 



0.6 0.8 

r/rO 



FIG. 3: Static quark potential at (3 = 1.3, 1/e = 1 on a 12^ lattice. The curve represents the 



continuum limit obtained by an interpolation of the results of 
with different orientations parallel to iJ's. 



28t |. Different symbols show y(r)'s 



IV. PERTURBATIVE RENORMALIZATION OF THE COUPLING 



In this section, we study whether the perturbative corrections are under control with the 
topology conserving gauge action. 

Two-loop corrections to the gauge coupling for general actions constructed by the pla- 
quette is available in j^^. Using that formula, the renormalized gauge couping qm defined 
in the so-called Manton scheme is given by 



— = — + Ai + A2g' 
9m ^ 



where the coefficients Ai, A2 are calculated as 
2N^ - 3 



S4- 



A2 = aR[sii2N^ - 3) + UN,{N^ + 1)] + 



15{N^ - 3N^ + 3) 



3(2iV2-3)(Ar2 + 3) 3,^. 



8N, 
2m 



— 2S4ti- 



{2 N^,-3){Nl + 2 



(10) 



-tl{Nl + l){Nl + 2). (11) 
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beta=1.3, 1/e=1 
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FIG. 4: Violation of rotational symmetry {V{r) — 14ont(|^))/K;ont(|^)! where Fcont(l^) denotes 
the continuum limit. Results for /? = 1.3, 1/e = 1 are shown. The error of V^ont(l^) is not taken 
into account (< 1%). 
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FIG. 5: {V{r) - V'cont(|r1))/V'cont(|^) for r = (1,0,0). Results with different f3 and 1/e values are 
plotted as a function of a/rQ. Left and right panels show the plot without and with the tree level 
improvement for the argument r, respectively. See the text for details. 
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Here, the parameters are A^^ = 3, S4 = -1/4!, sq = 1/6!, U = l/(4A^ce), te = l/{8N^e'^), 
uq = — l/(4!A'ce), and = —0.0011(2). Table IVl gives the next-to-leading and next-to- 
next-to-leading order coefficients Ai and A2 for various values of e. 

Since the perturbative expansion is poorly converging if one uses the bare lattice coupling, 
we also consider the mean field improvement using the measured value of the plaquette 



expectation value [31[ . To do so, we need a perturbative expectation value of the plaquette 
expectation value, which is available to the two-loop order for the general one-plaquette 
action j^^ as 



{W{1,1)) = 1-g' 



-W2{l,l)-g\N^-l)X{l,l) 



+9 



em 



-w,{i,ir-g'^^' 



QNr. 



-CZ(1,1) 



'12) 



Here the notations W^2(l; 1) and X{1, 1) are from the original calculation 



33I for the standard 



Wilson gauge action, and Z{1,1] 



1/V)W2{1, l)/4 (on a symmetric lattice V = L^) is 



introduced for generalization. Their values are ^^2(1, 1) = 1/8, X{1, 1) = —1.01 x 10 ^ and 
Z{1, 1) = 1/32 in the infinite volume limit. The constant C is written as 



C 



R 



d 



R 



(13) 



where C2(R) is the quadratic Casimir operator in a representation R of the group SU{Nc). 
dji denotes the dimension of the representation R, and T{R) is defined such that Ti jiitH^) = 
T[R)5°'^ for the group generator in the representation R. The coupling sji{P) is defined 
when we rewrite the gauge action in terms of a general form of the one-plaquette action 



1 



-^ReTi rP^Ax} 



(14) 



x,fi,u R 

where denotes the plaquette in the R representation. The values of these parameters 
for the topology conserving gauge action (jH) are 

11~ 

and T(3) = 1/2, T(6) = 5/2, T(8) = 3, C2(3) = 4/3, C2(6) = 10/3, C2(8) = 3. Using these 
numbers, we obtain C = 5 — 20/e and finally 

.2/5 5 



1 



(15) 



<^^^'^))-^-yni^-i44 



9 



(16) 
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Ai 


A2 


Bi 


B2 


Ai-Bi 


A2-B2 





-0.20833 


-0.03056 


-0.33333 


-0.03472 


0.12500 


0.00416 


2/3 


0.34722 


-0.04783 


0.11111 


-0.05015 


0.23611 


0.00233 


1 


0.62500 


-0.10276 


0.33333 


-0.13194 


0.29167 


0.02919 



TABLE V: Next-to-leading and next-to-next-to-leading order coefficients for the coupling renor- 
malization for various e. Mean field improved coefficients Ai — Bi, A2 — B2 are also shown. See 
the text for details. 

We define a boosted coupling as 

]l = L ( ^ : (l-^)A 

with the measured value of the plaquette expectation value P = {W{1, 1)) (see Table H)). It 
is defined to be a factor in front of F^^ when we rewrite P^^, = P explia^ F^^) and expand 
the action (jlj). The perturbative expansion of (|T7jl becomes 

k = -2+B, + B2g\ (18) 

g g 

where 

B,^J-(l-l\, + (19) 

3V ey V e/Vl8e 144; 9e 3e2 ^ ' 

We then obtain the perturbative expansion of the Manton scheme coupling in terms of the 
boosted coupling 

^ = 4 + (^1 - ^1) + (^2 - B2)f. (20) 

gii g 

Numerical values of Bi and Ai — Bi are listed in Table |3 We can confirm the effect of the 
mean field improvement; the two- loop coefficient A2 is significantly reduced by reorganizing 
the perturbative expansion as in ()20|). 

Using these results, the inverse squared renormalized coupling in the Manton scheme is 
obtained for each lattice parameters. In FigureElwe plot the coupling evaluated at a reference 
scale 5/ro as a function of lattice scaling squared. We use the two- loop renormalization 
equation for the evolution to the reference scale. Although the couplings are very different 
at the tree level, the one-loop results are already in good agreement among the different 
values of 1/e. Including the two-loop corrections, we find that the perturbative expansion 
converges very well and the agreement among different 1/e becomes even better. Good 
scaling toward the continuum limit can also be observed in this plot for the two-loop results. 
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FIG. 6: ^/ g\[ in the Manton scheme for various values of 1/e. The mean field improved expression 
(|20|) is used with the measured plaquette expectation value. Different symbols distinguish the value 
of 1/e (1 for circles, 2/3 for squares, for triangles). Open symbols with dot-dashed line represent 
tree-level results and open symbols with dashed lines are one-loop. The best results including 
two-loop corrections are shown by filled symbols with solid lines. 



V. STABILITY OF THE TOPOLOGICAL CHARGE 



In this section we discuss the stability of the topological charge with the topology con- 
serving gauge action. 

How the topological charge is preserved can be easily explained in the U(l) gauge theory 
in two- dimension, for which we can define an exact geometrical definition of the topological 
charge 



-E- 



27r 

X 

-iln(P^^(x)) 



(21) 



P^u{x) denotes the plaquette in the U(l) gauge theory. In two dimensions, Qgeo gives an 
integer on the lattices with the periodic boundary condition. The topological charge may 



change its value when the field strength pass through the point 
from F,',?* 



lat I 
fiu ' 



X] 



±7r. Since the jump 



X) 



-TT to Fl^fix] 



+71 is allowed with the usual compact and non-compact 



gauge actions, the topology change may occur without a big penalty. It is the U(l) version 
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of the Liischer's bound 

l-ReP^,(a;) <e (22) 

with e < 2, that can prevent these topology changes because the point Fji^^{x) = ±7r is not 
allowed under this condition. Furthermore, it can be shown that Qgeo is equivalent to the 
index of the overlap fermion with s = 0ife<l/5is satisfied. 

For the non-abelian gauge theories in higher dimensions, we do not have the exact geo- 
metrical definition of the topological charge (note that (jH)) gives non-integers). It is, however, 
quite natural to assume that a similar mechanism concerning the compactness of the link 
variables allows us to preserve the index of the overlap-Dirac operator for very small e. Also 
for larger e, we may expect that the topology stabilizes well in practical sampling of gauge 
configurations. 

Table. IVII summarized our data for the stability of the topological charge 

where r^i^^ is the autocorrelation time of the plaquette, measured using the method described 
in Appendix E of 3|. N^^^ denotes the total length of the HMC trajectories and is the 
number of topology changes during the trajectories. The topological charge Q is measured 
every 10-20 trajectories with the geometrical definition (jHl) after our cooling method. With 
this definition, Stabg represents a mean number of independent gauge configurations sam- 
pled staying a certain topological charge. But it only gives an upper limit, because the 
topology change is detected only every 10-20 trajectories and we may miss the change if Q 
changes its value and returns to the original value between two consecutive measurements. 
Therefore, our measurement of Stabg may give a good approximation when the topology 
change is a rare event. 

Results are plotted in Figure as a function of the lattice spacing squared. We find a 
clear trend that the stability increases for larger 1/e if the lattice spacing is the same. When 
the lattice size is increased from L/a = 12 to 16, the stability drops significantly for each 
value of 1/e. This is expected, because the topology change occurs through local dislocations 
of gauge field and its probability scales as the volume. For even larger volume {L/a = 20), 
our data are not precise enough, since the total length of trajectory is shorter. We also 
observe that the stability increases very rapidly toward the continuum limit. 
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Lattice size 1/e 


/? 


ro/a 




'^plaq 


#Q StabQ 






1.0 


3.257(30) 


18000 


2.91(33) 


696 


q 






1.2 


4.555(73) 


18000 


1.59(15) 


265 


43 






1.3 


5.140(50) 


18000 


1.091(70) 


69 


239 




2/3 


2.25 


3.498(24) 


18000 


5.35(79) 


673 


5 






2.4 


4.386(53) 


18000 


2.62(23) 


400 


17 






2.55 


5.433(72) 


18000 


2.86(33) 


123 


51 







5.8 


[3.668(12)] 


18205 


30.2(6.6) 


728 


1 






5.9 


[4.483(17)] 


27116 


13.2(1.5) 


761 


3 






6.0 


[5.368(22)] 


27188 


15.7(3.0) 


304 




16^ 


1 


1.3 


5.240(96) 


11600 


3.2(6) 


78 


46 






1.42 


6.240(89) 


5000 


2.6(4) 


2 


961 




9 /3 


2.55 


5.290(69) 


12000 


6.4(5) 


107 


18 






2.7 


6.559(76) 


14000 


3.1(3) 


u 


759 







6.0 


[5.368(22)] 


3500 


11.7(3.9) 


14 


21 






6.13 


[6.642(-)] 


5500 


12.4(3.3) 




90 




1 


1.3 




1240 


2.6(5) 


14 








1.42 





7000 


3.8(8) 


29 


64 




2/3 


2.55 




1240 


3.4(7) 


15 


24 






2.7 




7800 


3.5(6) 


20 


111 







6.0 




1600 


14.4(7.8) 


37 


3 






6.13 




1298 


9.3(2.8) 


4 


35 



TABLE VI: Stability of the topological charge Stabg. The length of the HMC trajectory N^,j, the 
autocorrelation time measured for plaquette Tpiaq, and the number of topology change are also 
summarized. 

For the study of the e-regime in a fixed topological sector, the lattices (1/e, /3, L) ~ 
(1, 1.42, 16) and (2/3, 2.7, 16) would be appropriate. Their physical size is L ~ 1.25 fm and 
the topological charge is stable for (100 — 1000)rpiaq trajectories. 
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FIG. 7: Stability of the topological charge. Different symbols represent the lattice size: L/a = 12 
(circles), 16 (squares), 20 (upward triangles). The value of 1/e is distinguished by the line type: 
1/e = (dotted), 2/3 (dashed), 1 (solid). Downward triangles are the data of [ll| measured on 
16^ lattices. 



VI. CONSTRUCTION OF THE OVERLAP-DIRAC OPERATOR 



A. Low-lying mode distribution of Hw 

We measure the low-lying eigenvalues of aHw on the gauge configurations generated 
with the topology conserving gauge action. We use the numerical package ARPACK js^, 
which implements the implicitly restarted Arnold! method. For the hermitian Wilson-Dirac 
operator aHy/ we take the form aHy/ = j^laD^ — 1 — s) with s = 0.6. 

Figure |H1 shows a typical comparison of the eigenvalue distribution for three values of 1/e 
on a 20^ lattice. The (3 value is chosen such that the Sommer scale ro/a is roughly equal 
to 5.3, which corresponds to a ~ 0.1 fm. From the plot we observe that the density of the 
low-lying modes is relatively small for larger values of 1/e. To quantify this statement we 
list the probability, P{< 0.1), to find the eigenvalue smaller than 0.1 in Table IVTTl For the 
above example, the probability is 41% for the standard Wilson gauge action (1/e = 0), but 
it decreases to 15% (9%) for 1/e = 2/3 (1). For another lattice spacing (ro/a ~ 6.5) and 
lattice size 16^, a similar trend can be found. In Table IVTTl we also summarize the ensemble 
average of the lowest eigenvalue Amin and the inverse of condition numbers Amax/Amin and 
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10 lowest eigenvalues of |H,„ 



10 



1e-05 



T 1 1 1 1 1 1 1 r 

Beta=1.3 ,1/e=1 + 
Beta=2.55 ,1/e=2/3 * 
Beta=6.00 , 1/e=0 □ 




_1_ 



_1_ 
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conf no. 

FIG. 8: Ten lowest eigenvalues of |a-ffvK| = {'y^^aDw — 1.6)| for gauge configurations with ro/a ~ 
5.3. Data are shown for 1/e = 1 (pluses), 2/3 (stars), and (squares). The lattice size is 20^. 

Amax/Aio, where Aio and Amax denote the 10th and the highest eigenvalues respectively. We 
may conclude that the lowest eigenvalue is higher in average for larger 1/e. 

B. Numerical cost 

In the numerical implementation of the overlap-Dirac operator one often subtracts the 
low-lying eigenmodes of aH^/ and treats them exactly. The rest of the modes are approxi- 
mated by some polynomial or rational functions. The numerical cost to operate the overlap- 
Dirac operator is dominated by the polynomial/rational part, because the subtraction have 
to be done only once for a given configuration. Here, we assume that 10 lowest eigenmodes 
are subtracted and compare the relative numerical cost on the gauge configurations with 
different values of 1/e. 

The accuracy of the Chebyshev polynomial approximation sgiiQ^^^^aHw) with a degree 
Npoiy can be expressed as \2^ 

{v\{l - sgnl^^^{aHw)f\v) 



{v\v) 



Aexp{-BNpo\y/K), 



(24) 



for a random noise vector |f ). A and B are constants. We find that they are A ~ 0.3 and 
B ~ 4.2 almost independent of the lattice parameters as shown in Figure El The reduced 
condition number k = Amax/ Aio enters in the formula with a combination A^poiy/ Therefore, 
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lattice size 


lA 


/3 


ro/a 


P{< 0.1) 


-^min 


-^min/ -^max 


-^lo/'^max 




1 


1.3 


5.240(96) 


0.090(14) 


0.0882(84) 


0.0148(14) 


0.03970(29) 




2/3 


2.55 


5.290(69) 


0.145(12) 


0.0604(53) 


0.0101(08) 


0.03651(27) 







6.0 


[5.368(22) 


0.414(29) 


0.0315(57) 


0.0059(34) 


0.02766(46) 




1 


1.42 


6.240(89) 


0.031(10) 


0.168(13) 


0.0282(21) 


0.04765(32) 




2/3 


2.7 


6.559(76) 


0.019(18) 


0.151(11) 


0.0251(19) 


0.04646(37) 







6.13 


[6.642(-)] 


0.084(14) 


0.0861(83) 


0.0126(15) 


0.03775(50) 


16^ 


1 


1.3 


5.240(96) 


0.053(13) 


0.111(12) 


0.0187(21) 


0.04455(31) 




2/3 


2.55 


5.290(69) 


0.067(13) 


0.1038(98) 


0.0174(16) 


0.04239(36) 







6.0 


[5.368(22) 


0.130(20) 


0.0692(90) 


0.0116(15) 


0.03451(62) 




1 


1.42 


6.240(89) 


0.007(5) 


0.219(13) 


0.0367(21) 


0.05233(26) 




2/3 


2.7 


6.559(76) 


0.020(8) 


0.191(12) 


0.0320(19) 


0.05117(29) 







6.13 


[6.642(-)] 


0.030(10) 


0.139(10) 


0.0232(17) 


0.04384(38) 



TABLE VII: The probability P{< 0.1) to find the eigenvalue lower than 0.1 for the hermitian 
Wilson-Dirac operator |ai?vy| = \j5{aDw — 1-6)|. Ensemble averages of the lowest eigenvalue and 
the inverse of condition numbers are also listed. The Sommer scale ro/a is the results for L = 16 
lattices. The values with [1 are from with an interpolation in /3. 

the numerical cost, which is proportional to A^poiy, depends linearly on k, if one wants to keep 
the accuracy for the sign function. From Table IVIII we observe that the reduced condition 
number is about a factor 1.2-1.4 smaller for 1/e = 1 than that for the standard Wilson 
gauge action. 

We also check that the above observation does not change by varying the value of s in 
a reasonable range. Figure ^1 shows a typical distribution of the low-lying eigenmodes for 
s = 0.2-0.7. We find that the advantage of the topology conserving gauge action does not 
change. Also, from these plots we can see that s ~ 0.6 is nearly optimal for all cases. 
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FIG. 9: The accuracy (^1)) as a function of Npoi/n for L = 16 (left) and L = 20 (right). We use 
4 gauge configurations and 10 values of A'poi = 60-195 for each parameter set. 
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FIG. 10: Ten near-zero eigenvalues for ro/a ~ 5.3 (left) and ro/a ~ 6.5 (right). Results are plotted 
as a function of 1 + s. Data for 1/e = 2/3 and are slightly shifted for clarity. 

C. Locality 

If the overlap-Dirac operator is local, the norm ||-D(a;, y)v(y)\ \ with a point source vector 
V at xq should decay exponentially as a function of x — xq jj] 



\\D{x,y)v 



C exp{—D\x — a;o|) 



(25) 



with constants C and D. This behavior is actually observed in Figure ^2 The plots are 
shown for different values of 1/e at the lattice scales ro/a ^ 5.3 (left) and 6.5 (right). We 
find no remarkable difference on the locality for different gauge actions. 

Recently, it has been pointed out that the mobility edge is the crucial quantity which 
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FIG. 11: {\\D{x,y)v{y)\\/\\D{0,y)v{y)\\)'^ with xq = measured on 10 gauge configurations for 
different values of 1/e. The lattice scale is rg/a ~ 5.3 (left) and 6.5 (right). 

governs the locality of the overlap-Dirac operator It would be interestina; 

to see the dependence of the mobility edge on the parameters in the topology conserving 
action, which is left for future works. 



VII. CONCLUSIONS 



We study the properties of the topology conserving gauge action (jlj) in the quenched 
approximation. For small e (< 1/20), the parameter to control the admissibility of the 
plaquette variable, it is theoretically known that the topology change is strictly prohibited, 
but we investigate the action with e ~ 0(1) for the use of practical purposes. With the 
(quenched) Hybrid Monte Carlo updation, we find that the topology change is strongly 
suppressed for 1/e = 2/3 and 1, compared to the standard Wilson gauge action. The 
topological charge becomes more stable for fine lattices, and it is possible to preserve the 
topological charge for 0(1,000) HMC trajectories at a ~ 0.08 fm and L ~ 1.3 fm. In the 
same parameter region, the standard Wilson gauge action changes the topological charge 
every 0(10) trajectories. The action is therefore proved to be useful to accumulate gauge 
configurations in a fixed topological sector. 

We measure the heavy quark potential with this gauge action at 1/e = 2/3 and 1. The 
lattice spacing is determined from the Sommer scale vq. With these measurements we also 
investigate the scaling violation for short and intermediate distances. The probe in the short 
distance is the violation of the rotational symmetry, and a ratio r^/ vq of two different scale 
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can be used for the intermediate distances. For both of these we find that the size of the 
scahng violation is comparable to the standard Wilson gauge action, which is consistent 
with the expectation that the term with 1/e introduces a difference at most O(a^). The 
action (@)) shows no disadvantage as far as the scaling is concerned. 

The perturbative expansion of the coupling and Wilson-loops is available in the literature 
for general one-plaquette action. We write down the coefficients for our particular action Q 
and observe that the convergence is very good if the mean field improvement is applied. The 
coupling constant in a certain scheme at a given scale is consistent among different values 
of 1/e. 

As a result of the (approximate) topology conservation, the low-lying eigenvalues of the 
Wilson-Dirac operator in the negative mass regime is suppressed. This is an advantage in 
the construction of the overlap-Dirac operator, since the numerical cost to evaluate the sign 
function is proportional to the inverse of the lowest eigenvalue for a given gauge configuration. 
In this case, the gain is about a factor 2-3 at the same lattice spacing compared to the 
standard Wilson gauge action. If the first several eigenmodes are subtracted and treated 
exactly, the gain is marginal, 20-40%. Similar improvements have been observed with the 
improved gauge actions, such as the Liischer-Weisz, Iwasaki and DBW2. 
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